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Abstract. We report on BeppoSAX simultaneous X- and y-ray observations of the bright y-ray burst (GRB) 990123. We 
present the broad-band spectrum of the prompt emission, including optical, X- and y-rays, confirming the suggestion that the 
emission mechanisms at low and high frequencies must have different physical origins. In the framework of the standard fireball 
model, we discuss the X-ray afterglow observed by the Narrow Field Instruments (NFIs) on board BeppoSAX and its hard X- 
ray emission up to 60 keV several hours after the burst, detected for about 20 ks by the Phoswich Detection System (PDS). 
Considering the 2-10 keV and optical light curves, the 0.1 - 60 keV spectrum during the 20 ks in which the PDS signal 
was present and the 8.46 GHz upper limits, we find that the multi-wavelength observations cannot be readily accommodated by 
basic afterglow models. While the temporal and spectral behavior of the optical afterglow is possibly explained by a synchrotron 
cooling frequency between the optical and the X-ray energy band during the NFIs observations, in X-rays this assumption only 
accounts for the slope of the 2-10 keV light curve, but not for the flatness of the 0.1 - 60 keV spectrum. Including the 
contribution of Inverse Compton (IC) scattering, we solve the problem of the flat X-ray spectrum and justify the hard X-ray 
emission; we also suggest that the lack of a significant detection of 15 - 60 keV emission in the following 75 ks and last 70 ks 
spectra, should be related to poorer statistics rather than to an important suppression of IC contribution. However, considering 
also the radio band data, we find the 8.46 GHz upper limits violated. On the other hand, leaving unchanged the emission 
mechanism requires modifying the hydrodynamics by invoking an ambient medium whose density rises rapidly with radius and 
by having the shock losing energy. Thus we are left with an open puzzle which requires further inspection. 

Key words, gamma rays: bursts - X-rays: bursts - radiation mechanisms: non-thermal 

1. Introduction 

GRB 990123 was one of the brightest y-ray bursts detected by the BeppoSAX satellite; it was also observed by the Burst and 
Trans ient Source Experiment (BATSE) on-board the Compton Gamma-Ray Observatory (CGRO) as trigger # 7343 jBriggs et alJ 
Il999t) . The Gamma-Ray Burst Monitor (GRBM) was triggered by GRB 990123 on 1999 January 23.40780 UT, approximately 
18 s after the CGRO trigger. The burst was simultaneously detected near the center of the field of view in Wide Field Camera 
(WFC) unit 1 IFeroci et alll999l) : the WFC "quick look" localized the burst wit hin a radius of 2 arcmi n (99% confidence level). 

With a redshift of z — 1.6004 and a luminosity distance of 3.7 x 10 cm feulkarni et alJll999al) . a y-ray (40 - 700 keV) 
flue nce of F v _ (1.9 ±0-2) x 10~ 4 ergs cm" 2 implies an isotropic energy release in y-rays alone of about 1.2 x 10 54 ergs (see 
also lBriggs et alJll999l:lKulkarni et alJll999al) . GRB 990123 would have been notable even just for this reason. Furthermore, it 
was the first burst from which simultaneous y-ray, X-ray and optical emission was detected. The prompt announcement of the 
burst position resulted in intensive multi-w avelen gth follow-up observations that brought a wealt h of new results: the discovery 
of prompt optical emission jAkerlof et al.ll 19991) . the detection of short lived radio emission jKulkarni et aljfl999bl) . the first 
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obser vation of a cle ar break in the optical afterglow light curve jKulkarni et alJll999allFruchter et alJll999HCastro-Tirado et alJ 
ll999llHolland et all2000l) . the first c onstraint on the polariz ation of a GRB afterglow JHiorfh et alJl999h . 

We refer the reader to the paper of Maiorano et al J fcooa for the complete analysis of the multi-wavelength afterglow data. In 
this paper we present the BeppoSAX GRBM and WFC observations of the y- and X-ray prompt event of GRB 990123 (section 
12. It and analyze the spectral properties of its X-ray afterglow (section l2~2l i. We discuss the observed prompt and afterglow emis- 
sion comparing them with the predictions of the "forward plus reverse shock" standard model. While the broad-band spectrum 
of the burst confirms a reverse shock origin of the optical flash (section l3~Tt . the multi- wavelength afterglow cannot be readily 
explained by basic models (section l3~2l . Errors on the parameters will be given at 90% confidence level (Ay 2 = 2.7 for a one 
parameter fit). 



2. The observations 

2.1. The prompt event 

The burst was detected by the BeppoSAX GRBM and WFC. The GRBM JAmati et all 997HFeroci et alJI 1997b consists of the 



4 anti-coincidence shields of the Phoswich Detection System, PDS ( Fronteraet alJl 997: Cos ta et all 19 98) and it operates in the 
40 - 700 keV energy b and. T he normal directions of two GRBM shields are co-aligned with the viewing direction of the WFCs. 
The WFCs tfageret all 19971) consist of two identical coded aperture cameras, each with a field of view of 40° x 40° full-width 
to zero response and an angular resolution of about 5'. The bandpass is 2 keV to 28 keV. The spectral resolution is approximately 
constant over the bandpass at 20% FWHM. In Fig. [flthe time profile of the burst is shown in various band-passes. The y-ray 
signal is detected for about 100 s with the GRBM (Fig. 1, bottom panel); it presents two major pulses, the brighter of which has 
a peak intensity of (1.7 ± 0.5) x 10~ 5 ergs cm -2 s _1 between 40 - 700 keV. The to tal fluence between 40 keV and 700 keV is 
(1.9 + 0.2) x 10~ 4 ergs crrT 2 ; this value improves the one given in lAmati et al.N2002l) . 

The burst profile at lower energies (2-10 keV and 10-25 keV) as measured by the WFC is shown in the top panels of 
Fig-E Near the end of the measurement the WFC was pointing close to the Earth horizon. Since the atmospheric absorption then 
plays an important role, the low-energy X-ray light curve was particularly affected: at 80 s after the trigger the Earth-atmospheric 
absorption is about 30% at 5 keV and the subsequent decay in the X-ray curve is partially due to the atmosphere. The portion 
of the light curve where more than 10% of the intensity is lost due to this effect is indicated by a dotted curve. The two X-ray 
light curves end when GRB 990123 sets below the Earth horizon. Note that the GRBM light curve is not influenced by the Earth 
atmosphere. The structure in the first 80 s of the X-ray lig ht curves indic ates a s oftening of each pulse, together with an increase 
in duration with decreasing energy. As remarked bv lBriggs et all ([1999) and by Frontera (2004), during the first two pulses the 
hardness of the GRB is correlated with the intensity. After the second pulse, a hard-to-soft evolution, typical of most GRBs, is 
seen. In the 2-25 keV and 2-10 keV energy bands, the fluence is (7.8 + 0.4) x 10~ 6 ergs cm -2 and (3.1 +0.3) x 10~ 6 ergs cm -2 , 
respectively; these estimates must be considered as a lower limits due to the eclipse discu ssed above. 

The model fitting of the data is done using the standard forward-folding technique (Briggsl l 19961) : in each instrument we 
assume a photon model and convolve it through a detector model to obtain a model count spectrum. The model count spectrum 
is compared with the observed count spectrum and the photon model parameters are optimized so to minimize a^ 2 statistic. The 
v f v data points are then calculated by scaling the observed count rate in a given channel by the ratio of the photon to count model 
rate for that channel; this ratio, and therefore the photon data points, are model dependent. 

The spectrum of the burst can be modeled by using the smoothly broken power-law proposed by B and et "all il 19931) . whose 
parameters are the low energy index 771, the break energy Eq, the high energy index 772. If we express E and Eq in keV, the Band 
function has the following expression: 

"Mi^f xexp (~l) (1) 

:exp( "HiooLvf (2) 



for E < (771 - 772) x Eq, and: 

\( m - m ) x£ () f-" 2 



N(E) . 



100 keV 



for E > (771 - 772) x Eq. In the case of GRB 990123, the best fit values for the 2 - 700 k eV integrated spectrum of the burst 
are: 771 = -0.89 ± 0.08, 772 = -2.45 ± 0.97 and E = 1828 + 84 keV JAmati et alJl2002l) . where E is expressed in the GRB 
cosmological rest-frame; at a redshift of z = 1 .6, it corresponds to an observed value of Eq = 703 ± 32 keV. 

Using the Band function as the model spectr um, we performed a spectral fitting of the WFC and GRBM data at the epochs 
of the detection of the prompt optical emission (Ake rlofet all 19991) . indicated in the lower panel of Fig. 1 . In Fig. |5]we plot the 
time resolved X- to y-ray spectra and the three ROTSE data points; in Table 1 we report the fluxes and the best fit spectral indices 
in the WFC and GRBM. 

Since the value of Eo derived from the 2 - 700 keV integrated spectrum is consistent with being above the energy range 
of the GRBM channels, we set the high energy spectral index 772 equal to its best fit value for the mean spectrum, 772 = -2.45 
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Seconds since GRBM trigger 

Fig.l. Prompt burst profile of GRB 990123 at X-ray (WFC: 2-10 keV and 10-25 keV) and at y-ray (GRBM: 40-700 keV) 
energies. Count rate (counts s ) is given as a function of time after the GRBM trigger, i.e. 1999, Jan 23, 09:47:14 UT. Since 
the X-ray burst was close to the Earth horizon, atmospheric absorption plays an important role. At 80 s after trigger the Earth- 
atmospheric absorption is about 30% at 5 keV and the subsequent decay is partially due to the atmosphere. The dotted line refers 
to the part where we loose more than 10% of the intensity due to this effect. The two X-ray light curves end when GRB 990123 
sets below the Earth horizon. Note that the y-ray light curve is not influenced by the Earth atmosphere. Typical error bars are 
given in the left part of the panels. Ind icated in the bottom panel are the time of the BATSE trigger and the prompt optical 
measurements (ROTSE, lAkerlof et alll999l) . 



( Ama tTet alJl2002|) : the values we find for Eq in the three time intervals are reported in Table 1. As can be seen in Fig. [2 the 
spectrum of the burst shows a hard to soft evolution; the spectral photon index for the first ROTSE exposure is 771 = -0.57 + 0.06, 
harder than the subsequent values of 771 = -0.95 + 0.09 and 771 = -1 .0 + 0.1, for the second and last time intervals, respectively. 

We find no evidence for an X-ray excess during the first ROTSE observation, as su ggested bvlBriggs et al.Nl999t) . In fact, our 
value for the low-energy spectral in dex is consistent wi th 771 = -0.63 + 0.02, found bv lBriggs et alJ (119991) analyzing the 10 keV 
to 10 MeV spectrum (see Fig. 3 of iBriggs et al.lll999i) : from the presence of an X-ray excess below 10 keV, we would expect 
to find 771 < -0.63 for v < 10 keV, in contrast with our results. Moreover, performing a spectral fit by using a simple power-law 
and considering only the WFC data at energies below 10 keV, yields a best fit photon index of -0.44 + 0.28, consistent with the 
spectral index we find for the overall spectrum. 
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Fig. 2. Simultaneous multi-wavelength spectra derived at three times during the burst (ROTSE V band, X-ray and y-ray, see also 
Table 1). The ROTSE data points have been connected with a line to the best fit model at 2 keV to guide the eye. The solid line 
connects the points relative to the first ROTSE observation, the dashed line those relative to the second ROTSE observation and 
the dashed-dotted line those relative to the third one. The corresponding times of the spectra are given in Table 1 in seconds after 
the GRBM trigger. The data relative to the first (third) time interval have been shifted up (down) of a factor 1000. 



Table 1. Fluxes an d photon indices of GRB 990123 prompt emission measured during the first three ROTSE exposures 
llAkerlofetalJl9 99). 



X-ray and y-ray data simultaneously with ROTSE 


start" 


end* 


m v 


Photon index 


E 


2-10 keV flux 


40-700 keV flux 


s 


s 




m 


keV 


10~ 8 ergs citT 2 s -1 


10~ 6 ergs cirr 2 s -1 


4.38 


9.38 


11.70 + 0.07 


-0.57 ± 0.06 




3.2 


7.3 


29.58 


34.58 


8.86 ± 0.02 


-0.95 ± 0.09 




5.7 


1.3 


54.87 


59.87 


9.97 ± 0.03 


-1.0 + 0.1 


250+160 

^ -70 


4.8 


1.2 



" Start time since GRBM trigger. 
* end time since GRBM trigger. 



2.2. The X-ray afterglow. 

The BeppoSAX follow-up observation lasted from 199 9 January 23.65 UT until January 24.75 UT. A previously unknown, 
bright X-ray source designated as 1SAX J 1525 .5 +4446 ( Heise et alJl999t) . was detected by the LECS and MECS units, at right 
ascension a = 15h25m31 s and de clination 6 = +44°46'.3 (equinox 2000), with an error-circle radius of 50". Within 22" 
of this position, Odewah n et alJ (fl999) detected a fading optical transient (OT) at magnitude R=18.2. During the first 10 min, 
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Table 2. Fit results for GRB 990123 X-ray afterglow, fix and fi 1 ^ are the spectral energy indices for the synchrotron and the IC 
component, respectively (see equation (0); the galactic hydrogen column density Nh is fixed to its value at the burst position, 
N H = 2.1 x 10 20 cm 2 ; is the hydrogen column density locally to the GRB site (in units of 10 22 cm 2 ); x 1 is reduced for 
the degrees of freedom dof. Case [1] refers to a single power-law model spectrum with variable spectral index; case [2] refers 
to the same model spectrum of case [1], but with a fixed spectral index fix = -1.23 (see Fig. [3}; case [3] refers to a single 
power-law model spectrum with the fit performed only on the LECS and MECS data points (PDS data excluded); case [4] is the 
synchrotron plus IC model spectrum, with the spectral index of the synchrotron component fixed at fix = - 1 -23 and the one of 
the IC component at fi 1 ^ = -0.73 (see Fig.[5]and Fig. [6). F ^ n and are given at 1 keV, in units of 10~ 3 photons cirT 2 s _1 keV -1 , 
in [1], [2], [3] and at 2 keV, in units of //Jy, in [4]. 
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the 2 - 10 keV photon rate of 1SAX J1525.5+4446 in the MECS is (0.16 + 0.02) count s" 1 , corresponding to (1.38 + 0.14) x 
10~ n erg cirT 2 s _1 . 

The 2-10 ke V flux steadily d ecreases with time following a power law decay F x (t) = F x (6 hr) x (f/6 hr) gx AtJy, with 
ax = — 1 .46 ±0.04 ( Maio rano et al.l2005l) . Here t is the time measured with respect to the y-ray event. We refer to lMaiorano et alJ 
(2005) for a complete description of the WFC, MECS and PDS data extraction and analysis relative to GRB 990123 light curves 
in the 2-10 keV and 15 - 28 keV energy bands. 

The observed fluence during the observations was 3.1 x 10~ 7 ergs cirT 2 (2 - 10 keV) and 7.5 x 10~ 7 ergs cirT 2 (2 - 60 keV). 
Extrapolating to long times after the trigger (about 1000 s) this implies an X-ray afterglow energy fraction of at least ~11% that 
of the prompt X-ray emission (2 - 25 keV). 

In the first 20 ks of the observation, between 1999 January 23.6498 UT and 1999 January 23.8813 UT, the total measured 
photon rate in the PDS is (0.20 + 0.05) count s _1 , equivalent to (1.9 + 0.5) x 10" ergs cm' 2 s" 1 ( 1 5 - 60 keV). The possible 
contamination by an X-ray source located in the MECS field has been analyzed by Maiora no et al.l (|2005). While its presence 
affects the PDS data in the second and last part of the observation, it is comparatively negligible during the first 20 ks: for the 
first time, an afterglow emission was detected at energies as high as 60 keV. In the following 75 ks and the last 70 ks of the 
observation, the high energy tail of the afterglow emission (15 - 60 keV energy band) was not detected. 

Since the first 20 ks of the observation are those in which we have the greater statistics and where the PDS signal is present, 
we focus our attention on the spectrum of the burst relative to this first time interval. We fit the LECS (0.1 - 2 keV), MECS 
(2-10 keV) and PDS (15-60 keV) data points during the first 20 ks of the observation, by using an absorbed power-law model 
spectrum where we set Nh = 2.1 X 10 20 cm 2 , for the galactic hydrogen column density at the burst positi on, and account for 
the eff ect of absorption locally to the GRB site by adding (see Table 2). The Nh was modeled using the lAnders & Ebiharal 
( 1982) relative abundances. To provide for flux inter-calibration when fitting simultaneously the LECS, MECS and PDS data, the 
LECS/MECS and PDS/MECS normalization ratios have been varied in the range 0.7 - 1 and 0.77-0.95, respectively dFiore et alJ 
1999). The best fit value for the spectral index is f3 x = -0.82 + 0.10 (see Table 2), consistent with the fix = -0.94 + 0.12 found 
by M aiorano et alJ (|2005) and obtained leaving unconstrained the range of values in which the normalizations ratio of the LECS 
and PDS can vary, and considering the LECS data between 0.6 keV and 4 keV (see Fig. 3 of lMaiorano et all 2005). 

Excluding the PDS data and fitting the only LECS and MECS points we get fix = -0.89 + 0.10, which is not significantly 
different from the previous result. 

The best fit values for the spectral indices relative to the next 75 ks and to the last 70 ks of the LECS and MECS observation 
are fix = -1.2 ± 0.2 and fix = -1.0 + 0.2, respectively. These two values are consistent, within the errors, with the one obtained 
during the first 20 ks, so they do not give any significant evidence of spectral evolution, in agreement with the results found by 
Maioran o et alJ J2005I) . 

3. Discussion 

3.1. The prompt event in the "forward plus reverse shock" standard model 

The comparison of the optical, X-ray and y-ray data presented in Fig. [5] shows that most of the energy is emitted in the y-rays. 
The extrapolation of the high energy time resolved spectra to optical frequencies falls at least 2 orders of magnitude below the 
simultaneous optical measurements, indicating the presence of an unobserved break between the optical and X-ray bands. This 
suggests different physical origins for the emission mech anism at low a nd high frequ encies, co nfirming the idea of a reverse 
shock origin for GRB 990123 optical flash dSari & PirarJl999aHGalama et all 19991 iBriggs et alll999t) . 
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Table 3. Closure relationships in the standard synchrotron fireball model and their corresponding values for /3 opt = -0.75 ± 0.07, 
B x = -0.82 ± 0.10, aw. = -1.1 + 0.03 and a x = -1.46 ± 0.04. 
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<v x <v c 




[a] 




3/2/Sop, 
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[b] a x - 3/2yS x = 


[c] a op , - a x = 


[d]B opt -B x = 




0.02 


±0.11 




-0.23 ±0.15 


0.360 ± 0.050° 


0.07 ±0.12 












< v c < v x 




[a] 


#opt — 


3/2y8 opt 


= 


[b] a x - 3/2/? x - 1/2 = 


[c]a opl - a x - 1/4 = 


[d]j3 opt -/? x - 1/2 = 




0.02 


±0.11 




-0.73 ±0.15* 


0.110 ±0.050 


-0.43 ±0.12 c 



" Not consistent with the expectations at the 7.2<r level. 
* Not consistent with the expectations at the 4.9tr level. 
c Not consistent with the expectations at the 3.6cr level. 



The current standard model for y-ray bursts and their afterglows invokes a "fireball" in relativistic expansion, probably within 

a collimated structure (jet). The high energy burst is thought to be produced by internal shocks developing from collisions of 

plasma shells with different velocities. When the outflow runs into the interstellar medium (ISM), it sweeps up the surrounding 

gas and heats it. A reverse shock does the same to the ejecta. In these shocks, the accelerated electrons radiate via synchrotron 

emission. The forward external shock, interacting with the ISM, produces the multi-wavelength afterglow. 

Since the shocked region behind the reverse shock is denser and co oler than the region behind the forward shock, it can 
j II -J | II I 

radiate in the optical band (Meszaros & Rees 1997; Sari & Piran 1999b). Optical emission from the reverse shock is expected 

to be detected during or soon after the high energy event dSari & Piranlll999al) . and was actually observed in GRB 990123 

llAkerlofetalJll999i) . Also the radio flare, observed at less than 3 d after the GRBM trigge r, can be interpreted as radiation 

from a reverse shock (Kulkar ni et JlT9 99b) rather than afterglow from a forward shock (Gala ma et all 19991). The BeppoSAX 

measurements of the prompt event agree with this scenario. 

3.2. GRB 990123 afterglow 

3.2.1. The optical to X-ray temporal and spectral indices in the standard fireball model 

In the standard picture dSariet all 1998). for a spherical shock of energy E propagating into a surrounding medium of density n, 
the afterglow emission F Y (t) scales as follows: for v above the cooling frequency v e , 

F y (t) cc £(p+2)/V/ (3) 
with a = -3p/4 + 2/4 and B - -p/2 and, for v < v c , 

F v (t) oc £^+WV /2 fV (4) 

with a = -3p/4 + 3/4 and B - — (p - l)/2. Here p is the power-law index of the shocked electrons. 

As we said in section l2~2l the afterglow light curve of GRB 990123 in the 2-10 keV energy band is well described by a 
power law of index ax - — 1.46 ± 0.04. During the first day, starting from 4. 1 hr after the burst, also the optical fl ux in the Gunn 
r band d ecreases steadily, and can be described as F, — 70 x (t/6 hr)"°p' u Jy with aw = — 1.10 ± 0.03, as found by |Kulkarni et alJ 
lll999al) : other autho rs give aw - -1.1 2 ± 0.08 (Holland et alJbOOOi) . or opt = -1.13 ± 0.02 dCastro-Tirado et al.lll999l) and 
a op t = -1.09 ± 0.05 dFruchter et al J 19991) . which are all consistent with the val ue of -1.10 + 0.03 tha t we adopt here. Two days 
after the burst, the optical flux declined more rapidly, with a opt — -1.65 + 0.06 dKulkarni et alll999al) . This steepening has been 
ascribed either to a transition of the fireball to a non-relativistic phase (Dai & Lu 1999) or to the signature of th e detection of a 
relativistic jet dRhoadslll99"9l ISari et alJll999t lMe~szaros & Ree"slll999t iHuang et alJl2000albld: Iwei & Lull2000l) . Since both are 
estimated to take place approximately 2 d after the GRB start time, i. e. after the end of the BeppoSAX follow-up observations, 
they would not be relevant for the X-ray afterglow, that can therefore be studied assuming a spherically symmetric relativistic 
expansion. Moreover, the faster decline at X-ray wavelengths is indicative of an evolution in a constant density environment and 
it is opposite to the expectations of the wind model dChevaher&Lill999l) . Thus, in the standard synchrotron fireball model, the 
optical and X-ray afterglow spectral and temporal indices of GRB 990123 should be compared with the closure relationships 
indicated in Table 3, that have been derived from equations and ©. 

We have seen in section l2~2l that the best fit value for the X-ray spectral index during the first 20 ks of the NFIs observation 
is Bx _ -0.82 + 0.10 . The power-law which best fits the optical-to-IR spectrum during the first two days has an index -0.8 + 0.1 
dKulkarni et al.lll999al) : this is consistent with th e spectral index measured from the optical spectrum of the transient 19 hrs 
after the burst: -0.6 9 ±0.1 dAndersen et al.lll999l) . Reconstructing the radio to X-ray afterglow spectrum on January 24.65 UT, 
Gala ma et ali i 1999) found a spectral index of -0.75 + 0.23 in the op t ical r ange and a spectral slope of -0.67 ± 0.02 between the 
optical and the X-ray wavebands; moreover, recently Maio rano et al.ld2005l) have estimated a best fit optical-to-IR spectral index 
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Fig. 3. GRB 990123 X-ray spectrum during the first 20 ks of the NFIs observation compared with the best fit power-law model 
of fixed spectral index Bx = -1-23. It is evident that the observed spectrum is too flat to be well reproduced by this model. 



of -0.60 + 0.04 . In the following d iscussion, we set y0 opt equal to the mean optical-IR observed sp ectral index of —0.75 ± 0.07, 
as derived by Holla nd et alJ (12000), and address the differences with the case f3 opt = -0.60 + 0.04 (Maiorano et al. 2005) where 
necessary. 

As shown in Table 3, if v opt < vx < v c , all the closure relationships are verified within 1 - 2 cr, except for the one relating the 
temporal indices of the optical and X-ray afterglow, that is not consistent with the expectations (~ 7<x level). If we consider that 
GRB 990123 is one of the brightest y-ray burst detected by the BeppoSAX satellite and that, consequently, the temporal indices 
are measured with relative errors less than 5%, this result enables us to exclude the case of a cooling frequency above the X-ray 
band. 

On the other hand, if v opt < v c < vx, the temporal slopes of GRB 990123 optical and X-ray light curves are readily explained 
within the standard synchrotron fireball model (Table 3). However, the closure relations involving Bx are not consistent with the 
expectations, giving evidence for a too flat X-ray spectral index. This last issue can also be addressed in the following way: the 
observed value of a opt allows to estimate that of p and, according to ®, gives p - -|ar pt + 1 = 2.47 + 0.04; from this we get: 
a x = QTopt - 1/4 = -1.35 ± 0.03,A, pt = |a opt = -0.73 ± 0.02, B x = fttopt - | = -1.23 ± 0.02. While the first two values are 
consistent with the observed ones within 2.2<x and lcr respectively, the latter is consistent with Bx = -0.82 ± 0.10 only at the 
4<x level. A spectral fit of the X-ray data points with a power-law of index fixed to Bx = -1.23 is shown in Fig.[3]and gives the 
results reported in Table 2: a similar spectral index is clearly not consistent with the data. 

It is then difficult to interpret GRB 990123 afterglow within the basic synchrotron model: while for the optical afterglow both 
the temporal and spectral behavior are well explained by simply assuming v c between the optical and the X-ray energy band at the 
beginning of the BeppoSAX. follow-up observation, in the X-rays this assumption only accounts for the slope of the 2-10 keV 
light curve, but not for the flatness of the 0.1 - 60 keV spectrum. 

Using the spectral indices found by Maio rano et al.l d2005l) . By — -0.94 + 0.12 and B opt = -0.60 + 0.04, the closure relations 
[a], [b], [c], [d] for this case are verified within 2.8cr, 3cr, 2.2cr and 1.3<x, respectively; thus, since the difference between the 
optical and the X-ray spectral index is enhanced, the observations are consistent with the relation [d], while [b] gives only 
marginal evidence for the existence of an X-ray excess. 
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Fig. 4. The U, B, V, R, and I band data at t = 1999 Jan. 24.497, 24.505, 24.496, 24.456, 24.487 UT JGalama et all 1999b. extrap - 
olated at t = 8.8 hr from the burst using a temporal decay index of o-opt = —1.1 and extinction corrected JGalama et alj fl999). 
are compared with the 0.1 - 60 keV spectrum corrected for galactic absorption (Nh = 2.1 X 10 20 cm~ 2 ) and for absorption 
locally to the GRB site (see the best fit res ults of case Til in Table 2). The dash- dotted, dashed and dotted lines are the power-law 
model s of fixed spectral index fi opt = -0.6 (Maio raiio et all 2005). A, pt = -0.75 JHolland et al.l2 000) 

> A) P t = _ 0.8 jKulkarni et alJ 

1999a), respectively, for the optical extrapolated data points; the solid line is the power-law model of case [1] in Table 2, where 
/3 X = -0.82. 



Finally, we can look at the optical to X-ray normalization. As already noted bv lKulkarnieTai] {l999a), at the epoch of the 
start of the NFIs observations the optical to X-ray spectral index appears to be rather fiat, y8x-opt = -0.54 + 0.02; this implies that 
the observed 2 keV flux is somewhat higher than one would exp ected from the extrap olation of the optical spectrum. As shown in 
Fig-El assuming a mean optical spectral index ofy8 opt = -0.75 (Holland et alJ l2000) the high energy extrapolation of the optical 
data falls well below the X-ray points; similarly, the low energy extrapolation of the X-ray data falls well above the optical ones , 
giving evidence for the presence of an X-ray excess. Assuming an optical spectral index of y6 opt = -0.6 jMaiorano et aljE005i) . 
the optical to X-ray normalization problem is solved by including a break at a frequ ency v = 1 keV a t t = 8.8 hr since the trigger. 
On the other hand, if one assumes an optical spectral index of /? opt = -0.8 (Kul karni et alJll999ah . the normalization problem 
is moreover strengthened, suggesting that the optical and X-ray emission could be related to different components. Different 
cases have been found for GRB afterglows in which there is evidence for an X -ray excess of t his kin d, and a contribution 
from IC scattering have been suggested as a possible explanation. In particular, lHarrison et alJ J 19991) analyzed the case of 
GRB 000926, finding that the broad-band light curves can be explained with reasonable physical parameters if the cooling is 
dominated by IC scattering; for this model, an excess due to IC appears above the best-fit synchrotron spectrum in the X-ray 
band. lCastro-Tirado et alJ (120031) analyzed the NIR/optical and X-ray spectra of GRB 030227, finding the two mismatching each 
other's extrapolation and suggesting that, in contrast to the NIR/optical band where synchrotron processes dominate, in the X-ray 
spectrum there could be an important contributio n of IC scattering. Finally, a mismatch in the optical to X-ray spectrum was also 
found for the bright GRB 010222 din'tZand et al.l200ll) . 
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Considering that in the standard model the values of the temporal indices are determined by the hydrodynamical evolution 
of the fireball, while the shape of the broad-band spectrum depends only on the assumption of synchrotron emission, from our 
analysis we conclude that: (i) using the spectral indices fix = -0.82 + 0.1 and /3 opS = -0.75 + 0.07, the optical afterglow is well 
explained by synchrotron emission and a standard hydrodynamical evolution (closure relation [a]), but the X-ray one seems to 
require an additio nal flat spectrum emission component (closure relations [b] and [d]); (ii) on the other hand, using the spectral 
indices found by iMaiorano et al.l J2005h . the shape of the broad-band afterglow spectrum is consistent with a single spectral 
component and standard synchrotron emission (closure relation [d]), but the temporal slopes of the optical and X-ray light curves 
are not consistent with the observations (closure relations [a] and [b]), suggesting the possibility of a different hydrodynamical 
evolution. We discuss the issue (i) in the following section and (ii) in section l3".2.4l 



3.2.2. A possible IC-scattering contribution: the synchrotron+IC model 

In order to interpret the combined spectral and temporal properties of GRB 990123 afterglow, we consider the possibility of 
relating the X-ray emission to an important contribution of IC scattering. 

The X-ray spectral index, fix = -0.82 + 0.10, similar to the optical one, ft opt = -0.75 + 0.07, and the difference between the 
optical to X-ray temporal slopes a opt - a% = 0.36, suggest the presence of an additional hard component. This naturally leads to 
a model in which the X-ray emission is dominated by IC-scattering of lower energy photons, while the optical one is synchrotron 
dominated; furthermore, an IC contribution could also explain the unusual high-energy (15-60 keV) emission observed by the 
PDS during the first 20 ks of the afterglow observations. 

In a synchrotron+IC model, the afterglow spectrum at some time to would be the sum of two power-laws: 



F v (/o) = exp [-ct(v)Nh - cr((l + z)v)N z H ] x 



UkeVi \2keW 



(5) 



if v < vj^(fo), where vjjf is the peak frequency of the IC 



with fix the spectral index for the synchrotron component; y6j, c = 
spectrum, otherwise, fi\^ = ^ ■ 

If the fireball expands in a medi um of con stant number density, if the expansion is spherically symmetric and if we are in the 
slow-cooling IC dominated phase dSari & EsintZOOll) . the contribution of IC scattering at a given frequency will evolve with time 
according to the following relations: 



..ic 



(6) 



Ffoct; 



1 a 



< V < 



(7) 



:r (9 P -H)/8 IC 



(8) 



plC K r (9p-10)/8+(p-2)/(4-p). 



(9) 



assuming a power-law approximation for the IC spectrum dSari & Esinl200ll) . 

Observing the above spectral and temporal relations, and considering the analysis we have done in section IT. 2. II it is clear 
that, if we want the synchrotron emission to be dominant in the optical band and the IC contribution to explain all the X-ray 
emission, we should have: (i) a cooling frequency of the synchrotron component above the optical band, since this condition 
explains well both the spectral and the temporal behavior of the optical afterglow (section l3.2.1> : (ii) a peak frequency of the IC 
component below the X-ray energy band, in order to account for the decreasing 2-10 keV light curve (if v)^ > 10 keV, the 
2-10 keV flux would increase with time, if dominated by IC emission). 

We can also notice that in the slow-cooling IC-dominated phase, the IC contribution affects the temporal slope of the syn- 
chrotron compone nt at frequencies above the cooling one, changing the temporal decay index from - 



3(p-l) _ 1 



to 



-3p z + 12p-4 
Mp-4) 



(see 



Sari & Esin 2001, equation (A. 20) in the electronic version of this paper). This means that the difference between the temporal 
decay index of the synchrotron flux at frequencies below and above the cooling one becomes j - instead of |: the presence 
of strong IC cooling makes the difference smaller than the case of pure synchrotron emission. Therefore, the case in which the IC 
contribution dominates the 2-10 keV emission seems to be favored with respect, for example, to one in which the IC component 
becomes dominant only in the high energy tail of the X-ray spectrum (say, for v > 10 keV), with the 2-10 keV emission still 



3(p-D 



-3p-+12p-4 



and 



being dominated by synchrotron emission. In fact, if this was the case, we should have a opt = - ~" 4 '' , ax - 
ttopj - axCopt + j^opt - jOx - j§ = 0. For GRB 990123, the latter has a value of 0.301 + 0.085, consistent with the expected one 



only within 3.5<x. 

For this reason, we analyze the hypothesis of the IC emission dominating the whole 2-60 keV spectrum. For v opt < v c < vx 
and vx > v 1 ^, with the synchrotron component dominating the optical emission and the IC component the X-ray one, the closure 
relations should be those indicated in Table 4, that for GRB 990123 are all verified within lcr. 
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Table 4. Closure relationships between the temporal and spectral indices in the synchrotron+IC model and their corresponding 
values forjSopt = -0.75 ± 0.07, /? x = -0.82 + 0.10, a opt = -1.1 ± 0.03 and a x = -1-46 + 0.04. 



v pt < Vc < vx and v x > vg 

a opt - 3/2j8op t = tt X - 9/4y3 x - 1/4 = a opt - 2/3a x + 1/6 = y e„ pl - J 8 x = 
0.02 + 0.11 0.13 ±0.23 0.040 ± 0.040 0.07 ±0.12 
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Fig. 5. GRB 990123 X-ray spectrum during the first 20 ks of the NFIs observation. The model spectrum is the sum of two power- 
law components (see equation 0): one represents the synchrotron contribution and has a fixed spectral index fix = -1.23, the 
other accounts for IC scattering and has a fixed spectral index f} 1 ^ = -0.73. 



Using the relations a opt = -1.10 + 0.03 = -3(p - l)/4 and a x = -1.46 + 0.04 = -(9p - ll)/8, we get for the index of 
the electron energy distribution p = 2.47 + 0.02 and 2.52 + 0.03, respectively. Using the spectral indices /3 opt = -0.75 + 0.07 = 
-(p - l)/2 and y8 x = -0.82 + 0.10 = -(p - l)/2 we obtain: p = 2.50 + 0.14 and p = 2.64 + 0.20, respectively. These values 
are all consistent with each other within ~ lcr, as expected from the closure relations being verified at the lcr level. We then 
choose the one with the smaller error, that is p — 2.47 + 0.02, derived from the optical temporal index aopt. In this way, we 
estimate the expected spectral indices between 0.1 - 60 keV for both the synchrotron component and the IC one; these are: 
y6x = 2/3a opt - 1/2 = -1.23 + 0.02, and/3£ = 2/3a p t = -0.73 + 0.02, respectively. 

In Fig.[5]and Fig.|6]we fit the 0.1 - 60 keV data during the first 20 ks of the NFIs observation with a model spectrum of the 
form (|5}, where we set fix = -1-23 and fi^r = -0.73. The best fit values (see Table 2) for the synchrotron and IC flux at 2 keV 
are F^ 1 " - 0.15*?'}? [i]y and - 0.72*j™ yuJy (x 2 - 35.5/26), respectively, consistent with the hypothesis of an IC component 
dominating on the synchrotron one in the X-ray band: F^/F^ < 0.15/0.72 = 0.21, taking into account that for the synchrotron 
flux at 2 keV we get only an upper limit. 

Our following step is to relate v)^, v c , f^kev anc ' at t — 0-37 d to the fundamental parameters of the fireball model: the 
fraction of the shock energy density that goes into the electrons and magnetic energy density, e e and eg, respectively; the energy 
of the fireball E52 in units of 10 52 ergs; the ambient medium number density ri\, in units of particles/cm 3 . 
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Fig. 6. GRB 990123 F v X-ray spectrum: the data and the model are those of Fig.|5]corrected for LECS/MECS and PDS/MECS 
flux inter-calibration normalization ratios and for absorption. 



Using the formulas given bv lGranot et ail Jl999h . lSari etail i ll 9981) . and bv lSari & Esii] l l200ll) in the slow-cooling IC domi- 
nated regime (see the Appendix in the electronic version of this paper for details), for the synchrotron cooling frequency we can 
write: 

v c (Hz) = A Vr (e B £ 52 )X< (10 ) 
where: 

A Vc (Hz) = (2.7 x 10 12 )(1 +z)5o£4)(170)^f^ 4 ' (11) 
with t4 the time since the burst measured in the observer's frame in units of days, and: 

"=20^4) <12> 

b " ih <13) 

For the peak frequency of the IC component: 

^(Hz)=A v ,ce^4^ 4 »- 1/4 (15) 
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where: 

A v ,c(Hz) = 2(2.9 x 10 15 )(1. 116 x 10 4 ) 2 (1 +z) 1/2 |^|j f rf 9/4 (16) 
For the peak flux of the IC component: 

F£Vy)=A F ic (£ 5 2»i) 5/4 4 /2 ( 17 > 
where: 

A F K(jj]y) = 0.39(2 x 1(T 7 )(1.7 x 10 4 )(1 + z)d^t l J 4 (18) 

with c/28 the luminosity distance of the source in units of 10 28 cm. For the flux of the synchrotron component at a frequency 

v > v c : 

FvWy) = A Fv 4n\E{ 2 4 (19) 
where: 

-1 I V — "2 2 -12 -3P 2+12 P~ 4 

A Fr (pJy) = (1.7 x 10 4 )(2.9x 10 15 ) £ !-(2.7 x 10 l2 ) l,2 v(Hz)- p/2 d^ I- — -j (170)^(1 + z)*F% 4< "- 4) (20) 
p 2 -2p-4 

d = ^T4T (21) 

(22) 

(23) 



2(p-4) 

f _ P 2 - 12 
7 4(p-4) 

(24) 

p-4 

Setting p = 2.47, z = 1.6004, d 2 8 = 3.7 dKulkarni et al]ll999al) . = 0.37 d, we have verified, by inverting the above relations 
(see the calculations in the Appendix), if reasonable values of eg, e e , £52, and ri\ can be found when v c > 2 eV and vj^(0.37 d) < 
2.0 keV. 

In order to minimize the energy requirement (see equation (A.61)), it is convenient to set (i) vj^(0.37 d) = 2 keV and 
(ii) v e = 2 eV, that is to say to set these frequencies at the higher and lower edge of the range in which they can vary to 
be consistent with the X-ray and optical data, respectively. These conditions are immediately evident taking into account that 
^Ikev = ^m'PkeV/vJJj') - ^, so F l £ can be expressed as a function of and F^kev me value of which is fixed from the 
observations; substituting into equation (A.61), yields £52 oc (v c ) flI3 (v^) au ~ ais , and for p — 2.47: £52 oc v f (v^)~ a36 . Thus, 
minimizing v c and maximizing gives the smallest value for £52. 

With the assumption (i) we have F)^ = F^ eV ; moreover, because of (ii), the synchrotron contribution at 2 keV extrapolated 
from the Gunn r data should be Ff" = 70(8.8 hr/6 hrT 1 46 (2 keV/2 eV)" 1 23 =s 10~ 2 /Jy, which is consistent with the X-ray 
data best fit value of £ ( ) Yn = 0.15^9 Jj yuJy within the errors. 



Fitting the 0. 1 - 60 keV spectrum with model Q where f3 x = - 1 .23, ft£ = -0.73 and £ ( , yn = 10~ 2 yuJy, yields Ff = Ff % 
F l £ = 0.83 ± 0.04 fjjy. We thus set v c = 2 eV, £^ yn = 10~ 2 /Jy, v% = 2 keV, F)£ = 0.83 //Jy, that give: 



2 keV 



e B = IQ- W e e = IQ- 1 (25) 

£ 52 = 7 x 10 4 m = 10 2 (26) 

With those values for the intrinsic parameters of the fireball, the temporal and spectral properties of the optical to X-ray afterglow 
are explained within the fireball model and the optical to X-ray normalization problem is solved. 

A similar calculation can be done setting £ ( , yn = 0.15 pJy and F^ = 0.72 /Jy, that are the values we get from the X-ray data. 
This choice yields a higher value for e# and a lower one for n\. 

e B = 6x IQ-* f ,s5x 10~ 2 (27) 



£52 = 6 x 10 4 rti = 8 



(28) 
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but has the problem that it predicts an optical extrapolated flux about 15 times above the observed one (0.15 /iJy/10~ 2 pJy = 15). 
The hypothesis of a "gray absorber" could explain an optical flux lower than expected from the X-ray data, without affecting the 
optical spectral slope, and alternative scenarios to ex plain simultan eously the optical and X-ray data of some GRBs by invoking a 
more gray extinction have recently been proposed by Stratt a et al J ( 120041) . for those GRBs for which some evidence of absorption 
locally to the GRB site was found. In our case, we do not have a strong evidence for local absorption (see Table 2), but we can 
neither exclude this hypothesis at all. 

We can also notice that with the values d25i . being e e len > 10, IC scat tering dominates the total cooling over the whole rela- 
tivistic stage of the afterglow evolution, as underlined bv lSari & Esin ( 2001 ); thus, the transition to the slow-cooling synchrotron 
dominated phase, in which the IC cooling rate is weaker than the synchrotron one and therefore has no effects on the synchrotron 
spectrum, should not be observable. This implies that for the 0.1 - 60 keV spectra relative to next 75 ks and last 70 ks, we do 
not expect changes in the spectral slope, in agreement with the results presented in section l2~2l However, we of course expect 
the flux level of the following spectra becoming dimmer of a factor determined by the slope of the X-ray light curve. Sinc e the 
15 - 60 keV upper limits are in fact consistent with a decay slope of a — -1.46 ± 0.04 (see Fig. 1 in lMaiorano et all2005l) . we 
thus conclude that the lack of a significant detection by the PDS in the following 75 ks and last 70 ks spectra, has to be related to 
poorer statistics rather than to a strong suppression of the IC contribution. 

Since we are before 2 d from the burst, we have used the isotropic dynamical equations and £52 is the fireball isotropic 
equivalent energy at the beginning of the afterglow phase. Considering that the observed isotropic equivalent y-ray energy, 
£iso(y) = 4-7iF y d 2 L ( \ + z) _1 , is 1.2 x 10 54 ergs, the efficiency rj y of the fireball in converting the energy in the ejecta E9 into y-rays 
has to be ^2 = r\ y s=s 0.2%, a reasonable value according to lPanaitescu et all i 19991) and Belobor odovl J2000h . 

To estim ate the energy of the jet £j at the beginning of the afterglow phase, we need the value of the jet opening angle 6 y 
According to lFrail et aTl l l200ll) : 



ft = 0.057 



(4T 



1 +z 



-3/8 



0.2 x £° 

ISO 



10 53 ergs 



-1/8 

(o.lcm- 3 ) 



(29) 



2.47 in our case), more than one power of t. However, the edge of the jet begi ns to be seen when y drops b elow the inverse 
spening angle l/0j. This occurs well before the sideways expansion starts t Panai tescu & Meszaroslfl999i) and the latter is 



Since fj = 2 d, we obtain #j = 0.064 = 3.7°; this implies £j = 10 54 ergs for the energy in the jet. 

This estimate on the j et opening angle is based on the simple relation y(t\) = ffr 1 , that itself cannot dis tinguish between 
the structure of GRBs jets llMeszaros et alJl998tlPai & Goul200lllRossi et alJ2002HZhang & Meszarosl2002l) . Conventionally, 
the late time temporal index 02 for th e optical af terglow of jetted GRBs, is believed to be the same as p. However, there are 
some caveats on this assumption (e.g. fWu et al.l20 04). most importantly the fact that the ambiguity of the un derstanding on the 
sidewa ys expansion of the jet leads to a great uncertainty on the value of a2- More specifically, as explained by Mesza ros & Reed 
l l 19991) for the case of GRB 990123, if one assumes th at the steepeni ng in GRB 990123 optical afterglow from f~ u to t~ h6S 
is caused by sideways expansion of the decelerating jet (Rhoads 1997), one expects a large steepening from t to t~ p (where 
P 

jet opening 

unimportant unti l the expansion is almost non-relativistic: GRB 990123, as well as other GRBs like GRB 010222, 020813, 
021004, 00091 1 dWu et all20 04). are good candidates for non-lateral expansion jets, from which we expect a steepening by t ' 
in the late time optical light curve. The observed difference of a\ - at = 0.55 + 0.07, matches the expectations at the 2.8<x 
level and is thus in agreement with the expected change in the decay slope from seeing the edge of the jet within a model of 
typical electron index of p — 2.5 (Meszaros & Rees Il999l) . These arguments are sufficient to justify the hypothesis of a jet 
str ucture within a s ynchrotron+IC model with p — 2.47, but there are of course several other plausible causes for steepening (e. 
g. iDai & Lull999l 

We can now compare our results with some previous analysis, bearing that those were not constrained by the 0.1 - 60 keV 
afterglow spectrum, which is analyze d in this work. For what conc erns the total energy in the fireball, a value of the order of 10 4 
for £52 has recently been proposed by Panaitescu & Kumar ( 2004) analyzing GRB 990123 early optical emission in the frame of 
the standar d "reverse-forward shoc k" model. 

In lp"anaitescu"& Kumarl fcOOlh . fitting the light curves in the radio, optical and X-ray energy band, it was proposed p = 
2.2H-2.4,e„ s9xlQ- 2 H-2xlQ- 1 ,e B = 10~ 4 H-3 x 10~ 3 , «, s 3 x 10~ 4 -=-2 x 10~ 3 , 6 S = l°-f3 and£j s 10 50 H- 5 x 10 50 ergs (values 



the isotro pic equivalent energy 
) found by Pana itescu & Kumarl 



taken from Fig. 2 in IPanaitescu & Kumarl200ll) : in particular, for £j = 5 x 10 50 ergs and Gj — 3 
is of ab out 4 x 10 53 ergs, 30 times less than the lower limit on the isotropic energy (E > 10 55 ergs) i 
(2004) and five times less than the observed isotropic energy emitted in y-rays. A value of 4 x 10 53 ergs requires the existence 
of a radiative evolution phase before the adiabatic one, causing the emission of more of the initial e xplosion energy an d leaving 
less for the adiabatic phase. The same considerations can also be done for the model proposed bv lWang et alJ J2000l) . who set 
p = 2.44, e e = 0.57, eg - 3.1 x 10~ 3 , n\ = 0.01, £52 = 5 and considered only synchrotron emission: in this case £52 is about 30 
times less than £i so (y). 

Finally, about the power-law index of the elec tron energy distribution, values of p = 2.44 and p = 2.3 were found by 
I Wang et all J2000l) and IPanaitescu & Kumarl J200ll) respectively, similar with the one of p — 2.47 we suggest in our analysis. 
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In this work, however, to fit the observed flat 0.1 - 60 keV spectrum and give an overall interpretation of the multi-wavelength 
afterglow, we are adding a dominant contribution of IC up-scattering to the synchrotron emission process. 

3.2.3. The radio observations: problems of the synchrotron+IC model 

GRB 990123 has been largely disc ussed within the sta ndard fireball theory and a reverse shock origin has been assign ed to both 
the optical flash and the radio flare dSari & Piranll999btlGalama et alll999llBriggs et alll999tlKulkarni et alll999bl). The light 
curves in the radio, optical and X-ray energy bands have been well reproduced ( Panaitescu & Kumar 2001 ; Kulkarnie^lXl 999b ) 



and different sets of values for the intrinsic parameters e e , eg £52 and n\ have been proposed (Pa naitescu & Kumari2 001. 2004 
I Wang et alJ200ol) . 

In this paper we are adding another piece of information to the puzzle: the spectrum of the optical to X-ray afterglow extended 
to the high energy tail of the 0.1 - 60 keV energy band. As we have underlined in the previous sections, GRB 990123 appears 
anomalous and difficult to interpret within a standard synchrotron model when combining this spectral information with the 
temporal properties of the optical and X-ray light curves. This is what has brought us to study the hypothesis of a synchrotron+IC 
model with values for the intrinsic parameters d25l and J26b . 

We want now to verify if such a model is consistent with the 8.46 GHz data. The radio afterglow was observed from t > 1 d 
after the burst (here we are not considering the data relative to the radio flare, to which we assign a reverse shock origin, 
iKulkarni et alJl999bl) . To estimate the value of the 8.4 6 GHz flux at about 1 d from the burst, we need to calculate the value of 
the synchrotron self-absorption frequency. According to Granot et "aTl fcOOlh . the self-absorption frequency v a is: 



v fl (Hz) = 0.247 x 4.24 X 10 9 (1 + zT 1 | ^— — ; \ — ~\ e^e^E^n 3 ' 3 (30) 



p + 2 \ 3/5 (p-D^ .-u/j^/s 
x 3p + 2) p-2 

Using the values ( I25l > and ( I26> we get v a = 15 GHz. 

Further, substituting those values in equations (A.l) and (A.2) of the Appendix, we get v m (l d) = 13 GHz = v a and f^ n = 
23 mJy. Thus we have: 

/8.46 GHZ\ 2 

F 8 .46GHz(l d) = 23 mJy lSQHz = 7 mJy (31) 

that is about a factor of 30 above the observations, that are lower than 260 /iJy llKulkarni et al.ll 999b). 

Another way to explain the problem is the following: over the first two days, the most abundant photons are between the 
radio and the optical regimes; if these photons are up-scattered by the Compton process, then the required Thompson optical 
depth is t = fx/fm = 0.83 pJy/F^ n . The observed radio/optical fluxes in the case of GRB 990123 are lower than ~ 260 pJy 
l lKulkarni et all 1 999b). that imply r > 0.83/260 = 3 x 10 3 . On the other hand, the optical depth has the following dependence 
on the parameters of the fireball model: 

t = cr T nR Mx 10 _7 (£/10 54 erg) 1/4 (n/lcm" 3 ) 3/4 (f/6hr) 1/4 (32) 

which has a value of ~ 10~ 4 , for £52 = 7 x 10 4 , ti\ = 100 and t = 8.8 hr. Thus, to reconcile this value of r with the observations, 
one should have F^ t n of the order of 3 x 10 _3 /10 4 x 260 pJy, that is to say 30 times higher than the observed one, that is the 
ultimate result of our synchrotron+IC model. 



3.2.4. An alternative solution for the puzzle 

In the previous sections we noticed that the observed temporal slopes of the optical and X-ray light curves of GRB 990123 
afterglow cannot be reconciled with the observed value of the X-ray spectral index during the first 20 ks of the NFIs observation, 
at least in the frame of the basic synchrotron standard model. Thus we tried to explain the X-ray spectrum by invoking a significant 
contribution of another emission mechanism in addition to the synchrotron one, that is IC scattering. However, we have some 
problems when comparing our results with the 8.46 GHz upper limits. 

A completely different approach could be to leave unchanged the emission mechanism, but change the hydrodynamics. 
Assuming that synchrotron radiation from electrons with a power law energy distribution is the only efficient mechanism, the 
shape of the afterglow spectrum is independent of the details of the hydrodynamics, thus the relation between p and fi is fixed. 
From the optical to X-ray spectral index, fix-opt = -0.54 + 0.02 dKulkarnietalll999al) . we can infer p — 2.08 ±0.04 if v c = 1 keV; 
thi s value of p predicts /3 n n t - -0.54 ± 0.02 and fix = -1.04 ± 0.02, that agree at the ~ lcr level with the spectral indices found 
by Maiora no et al.l ( 120051) . Assuming fi opt = -0.75 ± 0.07 and fix - -0.82 ± 0.10, the agreement is at the 2<j level. However, in 
this case, the closure relation yS op t -fix -1/2 = is only verified at the 3.6<x level. 

Further, in order to account for the temporal decay (a value of p — 2.08 predicts an X-ray temporal index of - 1 .06 ± 0.03 and 
an optical temporal index of -0.81 ± 0.03, that are not consistent with the observed ones), we need to change the hydrodynamics 
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and we do this by letting E oc f 5 * and n oc t 6 " where t is the time in the frame of the observer; we note that the standard model has 
8e — <5« = 0. Using and we find: 

a(v > v c ) = (p + 2)<W4 - 3p/4 + 2/4 (33) 
and 

a(v < v c ) = (p + 3)6 E /4 + 6J2 - 3p/4 + 3/4 (34) 

Substituting ax = -1.46 and ar opt = -1.1 into J33i and (I34> respectively, we obtain 6e = -0.39, so the shock should lose energy, 
and 6„ = 0.41. With those values for 6e and 6„ we obtain R(t) oc f 05 , where R is the radius of the shock. This result follows from 
the fact that R{t) ~ (Et/n) 1 ^ 4 . Thus, the density should increase rapidly with radius, as n oc R & . 

An increase in the density could let one think to the termin ation shock radius of the wind model bv lChevalier & Li| ([l999) or 
to the GRB jet running into a sharp-edged dense cloud. Indeed. lVietri e t al. ( 1999) suggest that the X-ray flare in GRB 970508 and 
GRB 970828 could arise from thermal bremsstrahlung emission of a shock-heated thic k torus surrou nding the central source. 
However, the observed X-ray spectrum is too steep to be compatible with this model. IShi & Gvukl i 19991) invoke a cloudlet 
to explain the 2-day radio flare (Kulk arni et al.l l 1999b). The energy demands of this model are of the order of 10 55 ergs, but 
ISari&Piranl (ll999b) have proposed a simpler model in which the radio flare arises in the same reverse shock which powers the 
prompt optical emission. Moreover, in addition to the necessity of justifying such a density profile, we also need to consider 
which observable effects could be expected from it e.g., the optical-near IR broad-band sp e ctrum of the burst becoming redder 
with time if the medium encountered by the fireball becomes very dense. Holla nd et all (|2000) have shown that there is no 
evidence for B opt varying with time between about 6 hr and 3 days since the burst, so that no observable reddening effect seems to 
become important during the afterglow observations. Using the relations n\ oc f - 41 and R oc f 005 , we can estimate that during the 
afterglow phase, while expanding for abouti?(3 d) -R(6 hr) = [(3 d/6 hr) 005 - l] R(6 hr) s 0.13 R(6 hr) s 10 16 cm in radius, the 
fireball should encounter a medium whose density increases of a factor (3 d/6 hr) 41 s 3; this means that starting with a typical 
value of «i = 1, the afterglow phase will entirely develop while the fireball is expanding in a region where the ISM number 
density is in the range of 1 - 3 particles/cm 3 . So, if we think the n oc R & profile only limited to the region of space where it is 
required to explain the afterglow observations, no extreme values of the ISM number density would necessarily be implied by 
this model, with no important reddening effects expected to be observed. 

The simplest way to lose energy is to assume that the shock is radiative. However, from the optical to X-ray normalization we 
know that v c = IkeV » v m and thus the shock is not radiative. A possible mechanism of ener gy loss was recently s uggested in 
relation with the issue of acco unting for bumps observed in the light curves of some GRBs (e.g. Scha efer et al.l2003tlKlose et al.l 
2004); SchaefeTHH] J2003I) proposed that the bumps observed in the optical transient of GRB 021004 could be related to 
inhomogeneities in the external gas in the form of clumps of denser material (Wan g et al J 120001) which should increase the 
afterglow brightness by enhancing the dissipation of the kinetic energy in the GRB remnant. This process could be of interest in 
constructing a physical framework for this alternative scenario, since we find the energy loss being connected with an increase in 

I I I L I I 

the density. Another possibility is energy loss via cosmic rays (e.g. Waxman 1995; Wick et al. 2004, and references therein). We 
know little about the energy carried off in cosmic rays in such strong shocks and so this possibility requires further inspection. 



4. Conclusions 

We have reported on BeppoSAX observations of GRB 990123 and discussed them in the frame of the standard fireball model. 
We analyzed the broad-band spectrum of the prompt emission, confirming the suggestion of a reverse shock origin for the optical 
flash. 

We have studied the properties of the 0.1 - 60 keV spectrum and compared the X-ray afterglow with the optical observations 
and the 8.46 GHz upper limits. The temporal slopes of the 2-10 keV and optical light curves are readily explained assuming 
v c between the optical and the X-ray energy band during the first 20 ks of the BeppoSAX follow-up observations, and setting 
p = 2.47 for the index of the electron energy distribution; this implies Bx = -1.23. Fitting the 0.1 - 60 keV spectrum with a 
single power-law model, yields Bx = -0.82 + 0.10, which is flatter than the value of Bx = -1.23 expected from the temporal 
slopes. 

We tried to relate the presence of this X-ray excess to the contribution of IC scattering and studied a model of synchrotron 
plus IC emission. By constraining ^kev ^m"' Vf anc ' v «? at tne ^ me °^ t ' le NFIs observations, we have found the corresponding 
values of the intrinsic parameters of the fireball model: with e e = 10 _1 , €b = 10~ 10 , £52 = 7 x 10 4 and ti\ = 10 2 , this model 
explains both the temporal and spectral behavior of the optical and X-ray afterglow, but violates the 8.46 GHz upper limits. 

We have compared and discussed our choice for the values of the intrinsic parameters, underlining the importance of the 
0.1 - 60 keV data as a new piece of information that should be taken into account when fitting the radio, optical and 2-10 keV 
light curves. 

Finally we have proposed an alternative scenario, where the problem of interpreting the observed value of both the spectral and 
temporal indices in the optical and X-ray energy band is solved by leaving unchanged the emission mechanism (only synchrotron 
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one) but modifying the hydrodynamics: in this framework, the totality of the observations could only be explained by invoking 
an ambient medium whose density rises rapidly with radius and by having the shock losing energy. 
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Appendix A: Details of the Synch rotron+IC Model 

The standard fireball theory for GRBs contains five parameters: the index p of the power-law electron energy distribution; the 
ratio e e of the energy in electrons to the post-shock energy in nucleons; the ratio e B of the magnetic field energy density to the 



post-shock nucleon energy density; the initial blast wave energy £52 (in units of 10 ergs) and the ambient number density n\ 
(in units of particles/cm 3 ). The value of the index p is generally estimated observing the decay slopes of the afterglow light 
curves or spectra. To find the values for the other four parameters, it is necessary to express them as functions of observable 
quantities, such as the ch aracteristic frequencies of the GRB spectrum or the amplitude of the observed flux at a given frequency. 
Wijers & Galama (1999) find the expressions for the intrinsic parameters of the fireball model as functions of the self-absorption 
frequency v„, of th e peak frequency v m , of the cooling frequency v c and of the peak flux F„ of the synchrotron spectrum. 
ISari & EsinI J200ll) find these expressions having also added the effect of IC scattering to synchrotron emission. 

In the case of GRB 990123 (see section 3.2 of this paper) we have estimates for v c , ^kev *»» anc ^ Thus, we need to 
express the intrinsic parameters e e , e#, £52 and n; as functions of those quantities. We start writing the direct relations that express 
Vc, F~ y , n .,, vJJr and F}!? as functions of e e , en, £52 and n\ and then we invert them. 

For the peak frequency and the peak flux of the synchrotron spectrum is (Gra not et all 1 999): 

v ffl (Hz) = (2.9 x 10 15 )(1 + z) 1 ' 2 l^j e l J 2 e 2 E l 5 ' 2 tf 2 (A-l) 

FT(piy) = (1.7 x 10 4 )(1 + z)4 /2 E 52 n\ /2 d 2 2 (A.2) 

where we have set pe akO?)/0peak(2.5) = 1 and tA T ,p a k(p)/iAp Pa k (2.5) = 1 in equation (25) and (26) of lGranot et alJ dl999t) . because 
those are slowly varying functions of p dGranot et al. Il999l) and we are interested to the case of p — 2.47. Note that dA.U is 



equivalent to equation (25) of lGranot et alJ dl999l) when setting f(p) = (^ry) and normalizing for /(2.5) i.e., 2.9 x 10 15 (1 + 

z) 1/2 /(2.5)(^) 2 ^V^V 2 = 3-22 X 10"(1 +Z ) 1/2 ^4 /2 ^V 2 ^no t et f ljl999|) 

In the case of slow cooling IC-dominated phase, for the synchrotron cooling frequency is (Sari et al. 1998; Sari & Esin 2001): 

v e (Hz) = (2.7 x 10 12 )e B 3/2 £ 52 1/2 «7'f rf 1/2 (l + tY^x' 2 (A.3) 

where x is the ratio of the the IC luminosity to the synchrotron one dSari&Esinl200ll) : 




(A.4) 

and d c is the time at which the slow cooling phase starts dSari & Esinl200i : 

t{f(d) = 170(1 + z)e 3 e B £ 52 «i (A.5) 

Note that (IA.3> is equivalent to equation (4.4) of ISari & Esii] d200ll) when expressing eg in units of 10~ 2 i.e., (2.7 x 

10 12 )(10- 2 r 3/2 (i^r 3/2 EJ^V 1 ^' /2 (1 + zr 1/ 2 x~ 2 _ 2.7 x 10 15 e ; V2 E^ 2 n- l t'/' 2 (l + zY^x' 2 iSari & Esinl200ll) . 
The peak frequency of the IC component is dSari & Esinl200 1): ' 

v£ = 2 T 2 v m (A.6) 
where y m is the minimum Lorentz factor of the electrons accelerated in the shock dSari & Esinl200ll) : 



y m = (1.116 xl0 4 )e^^-yj^j f/' a (A.7) 

Note that ( IA.7I > is equivalent to equation (4.6) in lSari & EsinI d200ll) when expressing e e in units of 0.5 and calculating (^ry) for 

p = 2.2, i.e. 1.1 16 x 10 4 x 0.5 fe) (j^f t~ vs = 93 e e 5 (— )' /8 tf' S dSari & Esinl200ll) . 

The peak flux of the IC component is dSari & Esin 200 1|): 

F^(pJy) = (2 X lO-^mR^FTiftJy) (A.8) 
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where fl 18 is the distance in units of 10 18 cm JSari et alJl998HSari & Esirl200ll) : 



^ 18=0 .39^) I/4 ? y 4 (A.9) 



In these formulas t& is the time measured since the burst in the observer's frame, in units of days, and J28 is the luminosity 

distance of the GRB source in units of 10 28 cm. 

Further, if we are in the slow cooling phase, the synchrotron flux at a frequency v > v c is JSari et alJl998l) : 



-(p-n 

,syn| Vc \ 2 ! r 



Fv = Fin — - (A.10) 
\Vtnl \Vc) 

We have thus collected all the quantities necessary to calculate the direct relations that express v c , F * yn , v 1 ^ and F)„ as 
functions of e e , eb, £52 and n\. Substituting \A.5\ in \AA\ and using the result in JA.3> we get the following expression for the 
synchrotron cooling frequency: 

v c (Hz) = A V( .(e B £ 52 ) fl «t< (A. 11) 
where: 

A„ c (Hz) = (2.7 x 10 12 )(1 +z)5^)(170)^f^ 4) (A.12) 

a = W^) (A - 13) 

6=—^— (A. 14) 
p-4 

c = l[ P —\ (A.15) 



Using jA.U and i n < IA.6I . we obtain for the peak frequency of the IC component: 
^(m^Ajc^El^ 4 (A.16) 
where: 

i5vi 1 v irvVn , ,u/2 (P~ 2 \\ -9/4 



A v ;c(Hz) = 2(2.9 x lO'^XLl x 10 4 ) z (l + z > \^ — [j f d (A.17) 
Substituting jA.9> and \A.2\ in JA.8> we get the expression for the peak flux of the IC component: 

F^QjJy) = A F ;c (£ 5 2»i) 5/4 4 /2 (A- 18) 
where: 

A F ic(/Jy) = 0.39(2 x 1(T 7 )(1.7 x 10 4 )(1 + z)c/ 28 2 ?y 4 (A.19) 



Finally, using equations ( IA. li . (IA.2> and (IA.3> in ( IA.10> . we get for the synchrotron flux at a frequency v > v c : 

F v (jiJy) = A Fv 4n[E f 52 4 (A.20) 
where: 

A F >Jy) = (1.7 x 10 4 )(2.7 x 10 I2 ) 1/2 (2.9 x 10 I5 )^v(Hz)^ /2 ,i 2 2 1 ^ — -j (170)^(1 + z)^>t d 4< "" 4) (A.21) 

= P 2 -2p-4 
4(p-4) 



P-2 
2(p-4) 



(A.23) 



p 2 - 12 

f = - (A.24) 
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p 2 -4p + 3 
P-4 



(A.25) 



The next step is to invert (IA.1 It . JA. 16b . JA. 1 81 . ( IA.20> and find the expressions for eg, e e , £52 and n\ as functions of p, v c , 
v~, F s J n and Fjp, the value of which we can constrain by using the observed spectral and temporal features of the burst afterglow. 
The results of our calculations are the following: 



e B = A^A^A^A^ (v c r (v^f (F^f (F v f 
where: 



a\ - 



a 2 = 



a 3 



I0p z - 


95/7 + 220 


5p 2 - 


■ 19p-4 


25(-p 2 


+ 5/7-4) 


2(5 P 2 - 


- 19/7-4) 


-5p 2 - 


31/7 + 204 



2(5p 2 - 19p-4) 

30(/7-4) 
5/7 2 - 19/7-4 



e e = A-«A-?A-;ZA-«« (v c y* {v^f (^ c f (F Y f 

where: 

/7 3 - 8/7 2 + 22/7 - 24 



a 5 = 



fl6 



a 7 



a 8 = 



2(5/7 2 - 19/7-4) 

4/7 2 - 17/7 + 4 
2(5/7 2 - 19/7-4) 

3/7-12 
2(5/7 2 - 19/7-4) 

-p 2 + 2/7 + 8 
5/7 2 - 19/7-4 



where: 

2/7 3 - 14/7 2 + 25/7 -4 



ag = 



aio = 



an = 



an = 



5/7 2 - 19/7-4 

/7 2 -5/7 +4 

2(5/7 2 - 19/7-4) 

5/7 2 - 17/7- 12 
2(5/7 2 - 19/7-4) 

_4p 2 +14/7 + 8 
5/7 2 - 19/7-4 



£ 52 = A^A^A^ A~y {v c f" {^) au {F^f (F v r 
-2/7 3 + IO/7 2 + 13/7 - 84 



an = 



a 14 



5/7 2 - 19/7-4 

9/7 2 -45/7 + 36 
2(5/7 2 - 19/7-4) 



(A.26) 

(A.27) 
(A.28) 
(A.29) 

(A.30) 
(A.31) 

(A.32) 

(A.33) 

(A.34) 

(A.35) 
(A.36) 

(A.37) 

(A.38) 

(A.39) 

(A.40) 
(A.41) 
(A.42) 

(A.43) 
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fll 5 = 



a 16 



5p 2 -p-76 
2(5p 2 - 19/7-4) 

4p 2 - 26p + 40 
5p 2 - 19/?-4 
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(A.44) 
(A.45) 



We now write the numerical expressions for (IA.26> . (IA.3H . ( IA.36> . JA.4H by substituting the values: z = 1.6004, c/28 = 3.7, 
p = 2.47, t d = 0.37 d, v = 2 keV: 



e^9.7xl0 5I (v e r (v^n^fOFvr 

where: 

fli = -2.27 

a 2 = -1.38 

03 = -2.37 

04 = 2.25 

e^2.2xlO- 9 (v,r(vL c f (F?)Vr)* 
where: 

a 5 = 8.32 x 10~ 2 
a 6 s 0.332 
a 7 = 0.112 
a 8 = -0.335 

where: 

a 9 = -0.121 

flio = 5.50x 10~ 2 
a u =0.575 
fli2 = -0.890 

E 52 - 1.1 x l<r 19 (v c r» (vL C f 4 (F^f 5 (F v )»« 

where: 

fli 3 = 1.03 

aw = 0.495 

d5 = 1.17 

flie = -8.99 x 10" 3 



(A.46) 

(A.47) 
(A.48) 
(A.49) 
(A.50) 
(A.51) 

(A.52) 
(A.53) 
(A.54) 
(A.55) 
(A.56) 

(A.57) 
(A.58) 
(A.59) 
(A.60) 
(A.61) 

(A.62) 
(A.63) 
(A.64) 
(A.65) 
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